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Abstract 

We show that in four-dimensional gauge theories dual to five-dimensional Einstein gravity cou- 
pled to a single scalar field in the bulk the derivative of the single heavy quark free energy in the 
deconfined phase is dFQ{T)/dT ~ — l/c^(T), where Cs{T) is the speed of sound. This general result 
provides a direct link between the softest point in the equation of state of strongly-coupled plasmas 
and the deconfinement phase transition described by the expectation value of the Polyakov loop. 
We give an explicit example of a gravity dual with black hole solutions that can reproduce the 
lattice results for the expectation value of the Polyakov loop and the thermodynamics of SU{3) 
Yang- Mills theory in the (non-perturbative) temperature range between Tc and 3Tc. 
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The conjectured duality between gauge fields and strings may hold the key to understand 



the complicated dynamics of non-Abelian gauge theories 



l|. The Anti-de Sitter/Conformal 



Field Theory (AdS/CFT) correspondence provides an explicit example where one 

of the main properties of A/" = 4 Supersymmetric Yang-Mills (SYM) in four dimensions, 
i.e., its exact conformal invariance, is already enough to determine that the background 
spacetime where the dual strings propagate contains a 5-dimensional AdS space. In this 
case, the dual string theory is a (critical, i.e., 10 dimensional) type IIB string defined in 
AdSs ® S5. Unfortunately, QCD and A/" = 4 SYM have very different properties even when 
Nc — > 00. For instance, while the former exhibits confinement below Tc and asymptotic 
freedom at high temperatures, the latter is always deconfined at any temperature in an 
infinite volume Moreover, in the region between and 4Tc lattice calculations have 
shown that violations from conformal invariance in QCD [5], or pure Yang-Mills (YM) with 
Nc > 3 [6j], are maximal. Thus, any realistic string theory dual description of QCD in 
the temperature range where strong-coupling effects are important must incorporate this 
important feature observed on the lattice. It remains unknown which type of critical string 
theory can describe the required infrared and ultraviolet properties of QCD, although some 
interesting top-down constructions have been proposed in the last years [7]. On the other 
hand, bottom-up phenomenological approaches have also been investigated in detail in the 
last years js, 9|. 

Significant progress can be made by considering 5-dimensional gravity theories in which 
the solutions of the equations of motion describe a spacetime that is asymptotically AdSs, 
which in turn implies conformal invariance in the UV for the gauge theory at the boundary. 
The minimal realization of this idea corresponds to Einstein gravity coupled to a nontrivial 
scalar field, ^ which is used to describe the violation of conformal invariance in the infrared 
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17|. In fact, it was shown in ll|] that there is a class of scalar potentials V{(f)) 



that generate black hole solutions with thermodynamic properties very similar to finite 
temperature QCD lattice results. According to the AdS/CFT duality dictionary, the non- 
Abelian 4-dimensional theories dual to these gravity theories are CFT's modified by the 
addition of a term ~ A^~^(9^, where A is the energy scale of the deformation and A < 4 
is the UV scaling dimension of the relevant operator dual to 0. One may interpret this 
gravitational model as a type of holographic "effective theory" for QCD valid at energy 
scales where effects from asymptotic freedom, which may not be fully described in the 
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supergravity approximation [3|, are negligible. Alternatively, such models appear quite 
naturally as pos sible realizations of non-critical string theories in 5 dimensions as discussed 
in detail in 10, isl- In a sense, currently gravity duals and perturbative QCD (pQCD) 
provide complementary descriptions of a plethora of QCD phenomena. 

Even though V{(j)) is not known from first principles, its form can be somewhat con- 
strained by requiring the dual gauge theory to display some important features observed in 
YM theories. In fact, there is a simple argument that justifies looking for potentials that 
have a simple exponential form V ~ —e'"^ [l^, where 6 is a positive number. The idea is 
to assume that the gravity dual's action comes from a non-critical string theory in D = 5 
dimensions and, thus, there should be a constant term that goes like a' (10 — -D) in the string 
frame action. This means that the Einstein frame potential V may have (at least) a term 



'act, the potentials studied in 



are such that V{(f)) ~ — e**^ in the 
13| that confinement occurs when (among other possible choices) 



that goes like e^'^l'^'^ . In 
IR and it was shown in 
V^(0) ~ — eV^"^ for large 0. 

Here we consider the model used by Gubser and Nellore in HI] and use the method 
(and conventions) they developed to solve the coupled nonlinear equations of motion for the 
metric, G^jy, and the scalar obtained from the action (in the Einstein frame) 

{d<\>f 



7^ 



-\/(0) 



(1) 



Conformal invariance in the UV can be obtained when — )■ at the boundary {u — > 0) 



and 



lim V{^) = + ^ A(A - 4)0^ + O(0^) 



(2) 



where the mass squared of the scalar m'^R'^ = A(A— 4) > —4 according to the Breitenlohner- 
Freedman (BP) bound [18]. Eq. (|2]) imphes that the spacetime will be asymptotically AdSs 
with radius R. 

In this effective theory, the gravitational constant k"^ is a free parameter expected to be 



symmetries is 



he most general ansatz for the line element compatible with the required 



= 0(m). 



(3) 



(note that has dimensions of length). The boundary of the bulk spacetime is located 
at u = 0. Finite temperature effects are included by considering solutions that display a 



horizon aX u = Uh where f{uh) = 0. The functions above are assumed to be finite at Uh- 
Using this ansatz for the metric the entropy density is given by 



. = 2vr^^ (4) 



while the temperature is (we use a prime for d/du) 



eAK)-^K)|f(^,)| 

^ = • 



We follow the gauge choice used in ll|] where (p = u. The equations of motion in this 
gauge are 

r-A'B' + l = 0, (6) 
o 

f" + i^A'-B')f' = 0, (7) 

GA'f + f{2AA'^ - 1) + 2e^^V = , (8) 

AA' - B' + ^ - =^V' = 0. (9) 

Gubser and Nellore showed that these equations can be solved using a generating function 
defined as G{u) = A'{u). In fact, 

A{u) = Ah+ f diG{i) (10) 

and, from ([6]) we obtain 



S,^^B,..j\<m±^ (U) 

and f{u) can be found from ([7]) and it reads 

f{u) = h rrfee-^^(«)+^(«) (12) 

where Ah, Bh, fi do not depend on u. Eq. ([H]) leads to an expression for V in terms of the 
quantities derived above 

-2B(«) / f'(v)\ 

V{u) = f{uf-^ [l - 24G{uY - 6G(n)^j . (13) 

The last equation of motion ([9]) does not add any new information at this point. 

This method to solve Einstein's equations becomes can be very useful because it is possible 
to show [11] that the equations above lead to the following nonlinear differential equation 
for G in terms of V 



Thus, the full solution of Einstein's equations for a given V can be directly obtained by 
solving the equation above for G. This can be done by choosing a value for Uh and then 
seeding a numerical integrator subroutine using a series solution of (HM around Uh- Such an 
expansion can be found using that Eqs. ([H]) and (P) when evaluated at the horizon become 

V{u,) = -3e-2^("'^)GK)/'(u,) (15) 
V'{u,) = e-2^(^'')/'(«,) , (16) 



which lead to the following identity 

nf^nA = - 



= -^TTTT-T (17) 



and the desired series 

Note also that Eq. f lT6|) implies that 

T = (19) 

This equation will be used in our discussion about the Polyakov loop later on. 

The main focus of this paper is to unveil some of the main properties displayed by 
Polyakov loops in deformed CFT's that have thermodynamic properties similar to 4 dimen- 
sional SU{Nc) pure YM theories. The relevant operator for our discussion is the path-ordered 
Polyakov loop defined as 

e{x) = — Tr P e'^o^^ M^,r)dT (20) 

is the non-Abelian gauge field operator and the trace is over the fundamental represen- 
tation of SU{Nc). Since in YM theories the only degrees of freedom are adjoint gluons, no 
other adjoint scalar fields (which would describe the transverse excitations of the set of Nc 
D3-branes) are included and we assume that the dual theory may be properly defined in a 
5 dimensional spacetime. 

When Nc ^ oo and the radius of the background spacetime is much larger than the string 
length, an infinitely massive excitation in the fundamental representation of SU{Nc) in the 
gauge theory is dual to a classical string in the bulk hanging down from a probe D-brane at 



infinity |19l-l2lj . The endpoints of tlie strings behave as non-dynamical fundamental probes 
and they are the analog of the infinitely heavy quark limit in quenched lattice gauge theory. 

According to the gauge/string duality, the correlator of Polyakov loops is given by 
(£(x)£*(0)) = Zstring, whcrc Zstring IS the full String generating functional that includes a sum 
over all the string worldsheets V (which may include connected and disconnected configu- 
rations) whose boundaries describe a QQ pair separated by a distance \x\ = L in the gauge 
theory j2l|. This defines the QQ heavy quark free energy via e~'^QQ^^''^^^'^ = {i{x)£*{0)). 
Within this approximation, the string dynamics (in Euclidean space) described by a given 
class of worldsheets V is given by the classical Nambu-Goto (NG) action 

SngCD) =tJ rfVVdet hf (21) 
Jv 

where Ts = is the string tension, hf = C^^d^Xt'd^X" {a,b = 1,2), X^^ = X'^(r,(j) 
is the embedding of the string in the 5-dimensional spacetime, and G^'^ is the background 
bulk metric in the string frame. In general, the metric in these frames are related as follows: 
= e^^C", where ^ > ^ 
The connected contribution to the correlator, {i{x)t{0)) - is ~ 0{N^) [21j]. In 

A/" = 4 SYM, the connected configuration that minimizes the action when LT ^ 1 is a 
U-shaped curve that connects the string endpoints at the boundary and has a minimum at 



some u* in the radial coordinate 
as (II]) was recently discussed in 



23| 



20[. The connected part of the correlator in models such 



In a confining theory, below Tc only connected worldsheets should be considered and, 
thus, in the supergravity limit the U-shaped configuration should be able to describe the 
full correlator at any L and T <Tc (other type of connected configurations become important 
above Tc at large L [21|]). In these theories, it can be shown that in the limit where L — t- cxo 
the vacuum heavy quark potential becomes Vqq(L) = ctqL + Eq + 0{1/L^), where 6 > 0, 
(To is the string tension, and the zero point energy Eq < [24] (the zero point energy in 



QCD was recently discussed in |25|). Roughly speaking, in the dual picture confinement 
occurs because the backreaction in the geometry given by the scalar prevents the bottom 
of the string to probe the deep IR region (as one pulls the quarks apart the only thing the 
dual string can do is to increase its length by the same amount). Thus, confinement implies 
that = when T < Tc. Similarly, the imaginary part of the Wilson loop induced by 



thermal fluctuations, which was recently derived in 



26| , can be shown to be identically zero 
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in confining theories when T <Tc. 

The disconnected piece of the correlator is the only 0{N^) contribution when Nc is large. 
In the gravity picture, this corresponds to two infinitely distant and, hence, independent 
straight strings going from the boundary of the space down to the horizon where the string 
endpoints can attach to the different D3-branes, which shows that this contribution 
indeed goes as N^. Note that the |x| — )■ oo limit of the full correlator is described by two 
different classes of worlsheets depending on T: below Tc one may use the connected U-shaped 
worldsheet (~ 0{N^)) while above cluster decomposition says that in this limit the full 
correlator equals the disconnected piece, which is ~ 0{N^). Therefore, in models such as 
confinement implies that jumps at the transition, which must be of first-order. Since 
the deconfinement phase transition in pure glue is expected to be first-order when Nc > 3, 
the class of gravity theories discussed here may indeed provide a description of the phase 
transition in YM theories at large N^. 

We shall now show how to describe the disconnected contribution to the correlator in 
theories such as ([I]). Using the gauge choice X^^ = (t, 0, 0, u), the regularized free 
energy of a single (static) heavy quark (up to a T independent constant) above Tc is 

Fq(T) =T,J^'du [^VM{u) - v/Mo(m)) - T, uf~' (22) 

where M{u) = \Gsoo\Gs^uu is a positive-defined function of u that is regular at the horizon. 
Near the boundary this function is Mq{u) = lim„^o^(^) = u ^-^ /(4 — A)^ as long as V 
obeys Eq. ([2]). The Polyakov loop is given by 

= e-FQ{T)/T_ (23) 



Instead of using Eq. ([22]), it is more convenient to work with its derivative dFq/dT 
{dFqlduh) duh/dT. In fact, 

dFc 



duh 

and one obtains that 



Ts^jM{uh) = Ts e^"" e^("'')+^("'') (24) 



dT Td\nT ^ ' 

The speed of sound in the plasma is = c? In T/c? Ins, which implies that duhjdXTvs = 

c^duh/dlnT. Therefore, one can rewrite the equation above explicitly in terms of 

dT ' Tel ds ^ ' 



On the other hand, one can use Eq. (jlj) to show that ds/duh = 'iA'{uh)s = 3G{uh)s. This 
identity, together with the horizon condition in Eq. f lT7|) and the formula for the temperature 
in (fT^ . shows that in gravity duals governed by Eq. ([T]) the following equation holds 

It - ^""^'ViMTDW) ' ^ ^ 

This is the main result of this paper. It shows that dFq / dT < (as expected from thermo- 
dynamics) and, consequently, that is a monotonically increasing function of T. For 
any potential that satisfies ^ one finds that in the small (high-temperature) limit 
Uh = {tcRT)^~^. Thus, in the high temperature limit Fq = —^T for any V (apart from a 
temperature independent constant and subleading terms that vanish at high T). This linear 
dependence in T is also observed in conformal theories 27|], which agrees with the fact that 
the spacetime considered here is asymptotically identical to AdS^. 

Since V{uh) is well defined at finite Uh-, one can see that dFg/dT only diverges when 
Cs — > 0. Moreover, for confining gravity duals the magnitude of the jump in \ {£) \ at (also 
known as the Gross- Witten point [28]) is related to how quickly vanishes at T T^. 
Continuous transitions where \ \ smoothly interpolates between zero and its high-T limit, 
which are expected to occur in the presence of dynamical fermions, can also be obtained 
depending on the choice for V{(f)). 

The class of potentials (which satisfies the BP bound in the UV and also leads to an 



interesting behavior at finite temperature) relevant for YM we use here is of the form 11 1 



= -12 coy + 6.^' ^ 

where 7 is directly related to the speed of sound in the IR and bi determines A near the UV 



fixed point. The choice 7 = y 1/2 and bi = 1.95 gives A = 3.3784 and leads to a black hole 
solution with the thermodynamic properties shown in Fig. 1, which can mimic the lattice 



data for the pressure p in SU{3) pure glue above as computed by Boyd et al. [29|] above 
Tc. This is done by computing the entropy density using Eq. (jl]) and choosing the value for 
/c| = 0.150376, which leads to the best agreement with the data. The other thermodynamic 
quantities are obtained from the entropy density via the usual thermodynamic identities. 

The gauge theory for this choice of parameters does not exhibit confinement (when 
^ = ^273) but has a very sharp (though continuous) transition near Tc, which is deter- 
mined by the minimum of c^. The trace anomaly 6 = 6 — 3p, where e is the energy density, 
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displays an interesting feature above Tc when divided by {TTcY: it plateaus in the temper- 
ature range between T/Tc = 1 — 4 where non-perturbative contributions to thermodynamic 
functions are expected to appear jsol . 31 1. At higher temperatures improved hard-thermal- 



32|. 



loop perturbation theory is expected to provide a good description of the data 

Because of our choice for V , the system is always in the deconfined phase and | (£) | can 
be computed at any T using ( 127|1 and there is no jump at Tc. The renormalized Polyakov 
OOP in quenched SU{?)) QCD was computed on the lattice in JbS] and more recently in 



34| . We compute |(£(T))| by integrating Eq. (!27|) and, because of the regularization, Fq is 



determined only up to a constant whose value is chosen via a match to the lattice data at 
T = 1.03 Tc (the first data point above Tc in Ref. 3J] for Nt = 4) for a given choice of a'. 
Our results for the renormalized Polyakov loop are shown in Fig. 2 (with a' /R^ = 3.5) and 
one can see they agree very well with the lattice calculations from [34] between T/Tc =1 — 3. 
At higher temperatures the curve starts to considerably deviates from the data, which is 
expected since the gauge theory considered here is not asymptotically free. 

In the class of dual theories considered here R'^/a''^ is not exactly the t'Hooft coupling as 
it occurs in = 4 SYM where the string dual is known. Here, its value could be fixed, for 
instance, by imposing that the ratio Tc/ y/d^ computed in a confining theory matches lattice 
results. However, for the non-confining theory discussed above R'^/a' is varied in order to 
find the optimal value that leads to the best description of the lattice data. Since this value 
turns out to be ~ C'(l), one should expect that higher-order derivative corrections to the 
gravity action, such as those studied in |35], are important. Because of these corrections, rj/s 
will acquire an explicit T dependence due to the coupling with (p evaluated at the horizon 
{sgI (without these extra terms rj/s = 1/ (Att) regardless of T IJiT])- In QCD, rj/s is expected 



to have a minimum near Tr 



38| . It would be interesting to find the minimal set of higher- 



order derivative corrections to ([T]) that can force this ratio to have a similar behavior in the 
dual gauge theory. 
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FIG. 1: Comparison between the pressure and the trace anomaly of the black hole solution of Eqs. 
([1]) and ()28p (solid lines) and the lattice data (dashed lines) from Ref. [2^. Only an interpolation 
of the lattice data is shown (the inclusion of the typical error bars would not lead to any significant 
changes in our analysis). 
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FIG. 2: The renormalized Polyakov loop (solid black line) for the gravity dual in Eqs. ([T]) and ([281 
and the lattice data for A'^t = 4 (red) and Nt = 8 (blue) from Ref. sj]. 
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